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I. Introduction 

A. Scope 

This report discusses the feasibility of making temperature profile 
measurements in the fuel preburner of the main engine of the space shuttle 
(SSME) using coherent anti-Stokes Raman spectroscopy (CARS). The princi- 
pal thrust of the work is to identify problems associated with making CARS 
measurements in high temperature gas phase hydrogen at very high pressures 
400 atmospheres). To this end a theoretical study was made of the char- 
acteristics of the CAR spectra of H 2 as a function of temperature and pres- 
sure and the accuracy with which temperatures can be extracted from this 
spectra. In addition the experimental problems associated with carrying 
out these measurements on a SSME at NSTLwere identified. A conceptual de- 
sign of a CARS system suitable for this work is included. 

Many of the results of the calculations made in this report are plot- 
ted as a function of temperature. In the course of presenting these 
results, it was necessary to decide whether the number density or the 
pressure should be treated as a fixed parameter. During a test 
firing, the SSME is throttled and the fuel preburner acts neither as a 
constant density nor a constant pressure environment. Since most physical 
quantities scale with number density rather than pressure (linewidth, op- 
tical dispersion, CARS signal strength, etc.) it was decided to present 
results holding density constant. The unit of density used throughout this 
report is the amagat (1 amagat density equals the density of a gas at 1 
atmosphere pressure at 0° C). Conversion from density to pressure is ac- 
complished via the equation of state. For convenience Appendix A presents 
isotherms for hydrogen gas for the pressures, densities and temperatures 
of interest in this study. 

B. Characteristics of the Fuel Preburner 

The following values were stated as design criteria for the tempera- 
ture profiling system in RFP No. 8-1-2-EC-25732 and are taken as the char- 
acteristics of the preburner where measurements are required. 

1. Temperature of gas: 

a. Steady state 2100 R° (1167 K) 

b. Transient maximum 3000 R° (1667 K) 

c. Normal minimum 1500 R° (833 K) 
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d. Transient minimum 200 R° (111 K) 

2. Gas composition from H 2 to Steam + H 2 + 0 2> 

3. Gas velocity: ^ 200 fps (61 m/sec) maximum. 

4. Chamber pressure: 5500 psi (374 atmospheres). 

5. Combustion chamber material: Haynes 188 + Zirconium Oxide 

Coating. 

6. Coating (wall) temperature: 2100 R° (1167 K) 

7. Measurement accuracy: 10 R° (5.6 K) 

8. Measurement resolution: 

a. Thermal - 10 R° (5.6 K) 

b. Spatial - from 0.25 to 0.50 inch square (1.6 to 3.2 cm 2 ) 

9. Time Response: 10 milliseconds 

Hydrogen and oxygen are fed into the fuel preburner in equal amounts 
by weight and hence lead to the reaction 16 H 2 + 0 2 -*■ 2H 2 0 + 14 H 2 . Pre- 
sumably this reaction takes place near the injector head and is completed 
by the time the gas mixture reaches the ends of the baffle plates and ports 
(see Figures 20a and 21b). Because little oxygen is expected to be present 
in the line of sight of the ports, it was decided at the initial contrac- 
tor's meeting that there would be little point in studying 0 2 as a probe 
species. In accordance with discussions held at Marshall Space Center and 
NSTL, the gas streamlines in the preburner are assumed to be steady except 
at the center. The maximum duration of a test firing is 500 seconds. In- 
formation on the fuel preburner, the SSME test stand, instrument accessi- 
bility to the preburner, physical environment adjacent to the preburner 
during a test run, etc. can be found in Section VII. 

C. Background 

CARS is a relatively new combustion diagnostic tool that is under- 
going rapid development [1,2,3]. It is an example of four-wave mixing in 
which three photons, usually two at a pump frequency and one at a 
Stokes frequency u) 2 mix through the third order susceptibility of a medium 
to produce a fourth or anti-Stokes photon at frequency u> 3 . The inten- 
sity of the anti -Stokes beam I 3 can be written 



where Ii and I 2 are the pump and Stokes laser beam intensities, z is the 
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length in which these beams coherently interact and x is the third order 
susceptibility of the medium which consists of a frequency dependent 
resonant part, x 1 + ix"> and a nonresonant part, x nr - When the frequency 
difference - u 2 approaches the frequency of a Raman active energy 
resonance of a major species in the bulk media, the resonant suscepti- 
bility becomes large compared to the non-resonant susceptibility. Because 
the intensity depends on the absolute value squared of the total suscepti- 
bility, the spectra appears as a resonant modulation on a non-resonant 
background. 

The resonant susceptibility is a function of the number density and 
the population difference between the levels involved in the transitions. 
This latter dependence means that the spectrum of a ro-vibration band of 
a molecule will be temperature sensitive since the population distribution 
of the rotational energy levels is governed by the Boltzman factor. Hence 
the temperature, number density, and bulk non-resonant susceptibility of 
the medium must be known in addition to the molecular parameters in order 
to theoretically model the CAR spectra of a molecular species in a multi- 
component gas medium. The bulk non-resonant susceptibility is obtained 
from a knowledge of the individual susceptibilities of the component gases 
and their approximate relative concentrations. The temperature can be ob- 
tained by examining the CAR spectrum of a major species. Comparison of the 
experimental spectra to calculated spectra has been the most reliable 
method of extracting accurate temperature information. Several excellent 
reviews of the theory and applications of CARS which discuss these consid- 
erations are available in the literature [4,5]. 

The accuracy with which temperature measurements can be made using 
CARS depends to a large extent on how well the CAR spectra of the chemical 
species under study can be theoretically modeled. The molecules present 
in the preburners are H 2 , 0 2> and H 2 0, all of which have been successfully 
modeled for atmospheric pressure CARS combustion studies [6,7,8]. How- 
ever, the complexity of the water spectra with its many closely spaced 
rotational lines invites one to consider this species for thermometry 
only as a last alternative. Since hydrogen will be the major species in 
the fuel preburner, it is the species of choice for temperature measure- 
ments. Oxygen is a trace species after combustion occurs and in accordance 
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with the above discussion will not be considered in this study. 

In many respects, is an ideal molecule for CARS thermometry. It 

£ _1 

has widely spaced rotational levels (B g = 60.85 cm ) reducing the num- 
ber of rotational energy levels populated at elevated temperatures and 
hence reducing the number of Q-branch spectral lines. In addition it has 
a large vibration-rotation interaction constant (a g = 3.06 cm"^) leading 
to widely spaced rotational lines. These properties mean that the in- 
dividual intensities of the few and widely spaced rotational lines can be 
studied which greatly simplifies analysis. Furthermore, the energy dif- 
ference between the lowest and first excited vibrational levels is so great 
(- 4400 crrf^) that no significant population of this excited level is ex- 
pected even at the maximum expected temperatures. The small linewidth 
pressure broadening coefficient (.0015 cm'Vamagat) and the large spacing 
between rotational lines (= 20 cnf^) indicates that spectral lines will 
not overlap even at the high densities present in the preburner. 

Thermometry using the broadband CAR spectra of H 2 has been performed 
in an atmospheric pressure H 2 ~air diffusion flame [6,7]. The relative 
intensity peaks of individual rotational lines were used to extract the 
temperatures in a range of 900 K to 2100 K. Standard deviations of 1 % - 8% 
were obtained with the poor precision due to the low signal-to-noise ratios 
that are characteristic of low number density environments. 
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1 1 . CcLtcJudLat&d CAR Specfrux 

There are several methods that can be used to extract temperature 
measurements from CAR spectra. Three are considered in this report. Two 
are simply variants of measuring intensity of individual ro-vibrational 
lines or groups of lines and these are discussed in Sections VI and VII. 

The third is to record all the rotational lines in a given spectral in- 
terval (i.e. 100 cm”^ ) and to compare this data with a library of com- 
puted spectra calculated for the expected experimental conditions and 
temperatures. This method is perhaps the surest and most accurate method 
of extracting temperature from the spectra since many lines are observed 
and anomalies can be easily detected. However, it is more difficult and 
expensive to implement than the other methods. 

Calculation of non-linear spectra, such as CARS, of a multiline 
molecular species is considerably more complex than the prediction of 
linear absorption or emission spectra. Much of this complexity is caused 
by the interference that takes place among the transitions and the inter- 
action of these transitions with the nonresonant background. A synopsis 
of the theory used to calculate CAR spectra is presented here. Various 
details can be found in the appendices. An effort has been made to point 
out some of the approximations used and limits of this theory. 

In general, the dipole moment per unit volume P induced in a medium 

-y. -y 

by electric fields E(r,i) can be written 
P(r,i) = i P (n) (r,i) 

n=l (1) 

= E x (n) Ui,a) 2 ,..o) n ) E(r,l ) E(?,2)...E(?,n) 
n=l n 

-y 

where r is the space coordinate, i represents the polarization of the field, 

. ( n \ 

u is the angular frequency and x is the susceptibility tensor of rank 
n + 1. If one considers only isotropic or centrosymmetric media, all even- 
numbered terms in this sum are zero because of inversion symmetry. Then, 
truncating the series at n = 3, Eq. (1) can be written as 

P a (r,1) = x^^wjE^Wi) + x^g| 5 ( w i ,!i> 2 >^ 3 ) E e (r,l)E^(r,2)E 6 (r,3) (2) 

where the summation over repeated indices convention is used. 
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In general, has 81 elements. However, for isotropic media only 
21 are nonzero, of which only 3 are independent [9] 


x nn 

X xxxx 

= x yyyy 

x zzzz 

x i 1 22 

x xxyy 

“ x yyxx 

x xxzz 

x 1212 

x xyxy 

" x yxyx 

x xzxz 

X 1 221 

x xyyx 

X yxxy 

x xzzx 

x llll 

= X 1 122 

+ X 1 212 

+ X 1 221 


x yyzz 


zzyy 


x yzzy x zyyz 


Most practical applications of CARS use the same frequencies for two 
of the E fields (three-color CARS) making x-|-] 2 2 = x 1212’ since x is 
invariant to the six permutations of the pairs (e,u>i), (y,u 2 ), (6,w 3 )^ 
Hence, for this case there are only two independent elements (i.e. x -|in 
y (3) ) 

x 1221 ' * / o \ 

If the summation of the x v ' of Eq. (2) is performed for these con- 
ditions, the third order polarization vector becomes [10] 

= jfi X^]22^“ 1,a i ) 

+ Sx^^-] (“l E(l) 2 E(2)* (4) 

x exp [i(2k 1 -k 2 )*r] 

- *CARS ®CARS E <’)W expCK^-^J-fl 

where the input fields E(l) and E(2) are assumed to be plane waves with 
polarization vectors e x and e 2 and wave vectors in the medium of ki and 

k 2 (k = nu/c, n = y/\ + 4ir X ^). If the Stokes and pump fields are par- 
allel (e 2 = ei), e CARS = e x (the three fields have the same polarization) 

X CARS = 6x 1122 + 3x l 221 = 3x llll (“1 » M l *““ 2 ) ( 5 ) 

If the polarizations of the Stokes and pump fields are perpendicular 
(ei_i_e 2 ), the anti-Stokes polarization is parallel to the Stokes field 

(e CARS “ e > Wlth c o v (3) (6) 


' CARS x cflRS ■ 3x,^ <-I. 

In this latter configuration, the CARS susceptibility is less than in the 
first configuration by the term 6xi||l* 
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Substitution of the linear polarization into Maxwell's equations 
leads to the usual plane wave solution with frequency and polarization 

vector unchanged but with wave vector k = where n = -y/l + 4ux^ ^ . Sim- 
ilarly, one can obtain an expression for the intensity of the CARS wave 
with Maxwell's equations using Eq. (4) for P^). If, in addition to the 
above restrictions, it is assumed that 

1. all input fields are plane monochromatic waves, 

2. the polarization vector is written as an expansion of plane 
wave frequencies, 

3. the CARS wave is a plane wave whose amplitude can be a 
function of longitudinal position, and 

4. Ak = 2k].-k 2 -k 3 <<2k 3 , 

then the time-averaged intensity of the anti -Stokes wave created in a length 
l is given by [10] 


= 16 u 4 gjg 

n 1 n 2 n 3 c l + 


lx 


CARS 


III 2^ 


rsin(Ak£/2) "1 2 

L Ak£/2 J . 


For exact phase matching Ak = 0 and 


(7) 


I 3 


1 6tT 4 0) 2 

n!n 2 n 3 c 4 



I?I 2 * 2 


( 8 ) 


The frequency dependence of x^rs is obtained from a quantum mechani- 
cal calculation of the microscopic molecular dipole moment arising from the 
third power of the microscopic fields. Averaging this over all molecular 

orientations and multiplying the result by the specie number density yields 

(3l (3l 

the bulk dipole moment P . x is obtained as the coefficient of the 
various third order products. Using diagrammatic perturbation theory one 
obtains 24, four-fold sums over all molecular states [11]. 

If 


1. all terms containing no Raman resonant terms in the denominator 
are lumped into a nonresonant term, v , 

NK 

2. there are no one- or two-photon resonances, and 

3. the Stokes and pump fields have parallel polarizations 
(ei = e 2 ), 

then Xpadc can be written in terms of the differential cross section for 
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spontaneous Raman scattering, for two stationary states (k and n) of 
the molecule, viz [10] 


+ !VVhi d g, i , 9 ) 

X CARS " x nr f\n 2 <4 W \ n -“ i + “ 2 - iY kn 

i, L 

where N. is the number density in the i state, is the energy dif- 
ference between states k and n, and Y kn is the homogeneous line half 
width (HWHM) . (Equations (7) and (8) differ by a factor of 1/16 and 
Eq. (9) differs by a factor of 2 from those given by Ref. [10].) 

As noted above, the equations were derived under the assumption that 
the input radiation is monochromatic. However, all real light sources 
have finite line widths, and in particular, CARS employed for thermom- 
etry on practical combustion systems uses a broadband Stokes field in order 
to simultaneously view many transitions of molecular species. Yuratich 
[12] has considered the problem of finite line widths of the incident 
fields in detail. If the input waves are approximated as stationary 
stochastic processes, he finds that the anti-Stokes spectral intensity 
I3 is given by the double convolution 


I3 ^i(^i)[/|x CA Rs(— ^ 


1-1 ( u 3 + i 02 -w l ) 1-2 (^ 2 ) *^2 


dwi (10) 


where Lx (»].) is the normalized lineshape of the pump intensity (using a 
single pump laser). L2 (^2. ) ^ the normalized lineshape of the Stokes in- 
tensity. 


Equation (9) is the expression for the susceptibility for a single 
resonance. A molecule has many closely spaced resonances and the CARS 
susceptibility will be the sum over all such resonances, i.e. 



(ID 


where N is the number density of the Raman active molecule and is the 
population difference between upper and lower vibration-rotation states. 
Substituting Eq. (11) into Eq. (10) and performing the integrations over 
frequency intervals where the molecular resonances occur will yield a CAR 
spectrum for realistic physical situations. 
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If the line width of the Stokes radiation is assumed to be much greater 
than the pump laser line width, Aw 2 »Aw 1 , then the narrow do^ integration 
can be performed and 


iU 3 H j x cars (w2-“i) 


Iq(uq)a. / |*,, nnc ((jj 9 -w?) 1 2 1 ^ 0 ) 3 + 012 -^ 1 ) L 2 (u) 2 )du 2 (12) 

where o>f is the center frequency of L^Wi). Using Gaussian lineshapes for 
the laser spectral profiles, this integral can be numerically performed to 
obtain a CAR spectrum. This is the approach used successfully by Hall [14] 
to generate synthetic CAR spectra of Ng and CO. 

It is possible to integrate Eq. (10) using Lorentzian or Gaussian line- 
shapes and obtain an analytical solution in terms of the complex error func- 
tion without resorting to the broad line width approximation. These solu- 
tions were obtained for several combinations of Gaussian and Lorentzian 
lineshapes by Yuratich [12]. 

Spectra presented in this report were calculated using a computer 
program based on the analytical solution of Yuratich. Comparisons of 
CAR spectra of molecular nitrogen using the two methods yield similar re- 
sults [13]. 

Evaluation of Eq. (10) begins with the calculation of xq^rs* Re- 
writing Eq. (11) results in 

X CARS 


x nr + ? 

J 


a 3 


Q .+u 2 -wi-iY i 

J vJ 


(13) 


where 


a j 


_ Nni c 4 

hn 2 uf" 


A 

A j dn 


The rotational selection rules for CARS are the same as those for 
Raman scattering [10], i.e. AJ = 0, + 2. For forward scattering, the 
spontaneous Raman cross section can be expressed in terms of the compo- 
nents of the derived polarizability tensor. For the Q branch [14] 


da _ 4 |i a 2 , _7_ 

dn ' c ' 2u ) 0 [ M 45 


j£b J l 

M D J 


(v + 1 ) 


(14) 


where M and u 0 are the reduced mass and molecular oscillator frequency 
respectively; bj is the Placzek-Teller coefficient [5]; a and y are the 
derivatives with respect to internuclear coordinate of the mean molecular 
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polarizability and anisotropy; v is the vibrational quantum number of the 

initial level. Evaluation of Eq. (14) for the case of H 2 is delineated 

in Appendix B. The calculation of the population difference A. of Eq. (13) 

J 

is outlined in Appendix C and the formulae for the energy differences ft. 

3 

are presented in Appendix D. The homogeneous line half widths, y., are 
functions of temperature and pressure and a discussion of these quantities 
is presented in Section III. x NR is a measured constant that must be 
scaled with number density. The value for the nonresonant susceptibility 
for deuterium was used in this study [17]. 

The lineshapes L(x) in Eqs. (10) and (12) may be of any profile. 

Laser line widths may be Lorentzian 



L L (x) " ir(x z +y ) z ’ 

Li 

f - L (x)dx = 1 

(15) 

or Gaussian 

exp(-x 2 /y 2 ) 

L G (x) ~ 

y g ^ 

; y*L G (x)dx = 1 

(16) 


where y_ = y/ yJln'L and y are HWHM. For multimode laser fields, the out- 
put is more closely approximated by a Gaussian profile. Substitution of 
Eq. (16) into Eq. (12) for the Stokes and pump laser fields and integrating 
yields broadband CAR spectra. This calculation using the analytical 
expression of Yuratich is detailed in Appendix E. 

Hg CAR spectra were calculated for a range of temperatures, densities, 
and Stokes beam widths. Figures (1) and (2) show computed, normalized 
spectra for two temperatures (800, 1600 K) that span the operating range. 
The four spectra in each figure show the effects of the variation of num- 
ber density from 1 amagat (upper row) to 100 amagat (lower row) and the 
variation of Stokes beam width from 100 cm -1 (left column) to 250 cm -1 
(right column). For all spectra, a Stokes beam center of 14679 cm" 1 was 
used, so that the resultant anti-Stokes beam is centered between the Q(0) 
and Q(5) lines. Appendix F contains a series of such spectra for 
800 < T < 1800 K, at 100 K steps. 

The prominent features of the spectra in these figures are the line 
broadening with increasing temperature and density, and the relative 
peak intensity rescaling with Stokes beam narrowing. Since the spectral 
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INTENSITY (NORMALIZED) INTENSITY (NORMALIZED) 




WAVENUMBER (cm" 1 ) 




Figure 1. Calculated CAR spectra of H 2 for a temperature of 800 K. Spectra were computed for 
densities of 1 and 100 amagat and Stokes beams with spectral widths of 100 and 250 cm" . Each 
spectrum is normalized to the maximum peak intensity for that spectrum. 
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Figure 2. Calculated CAR spectra of H 2 for a temperature of 1600 K. Other details same as those 
in Figure 1 . 


extent, Q(0) to Q(5), is about 100 cm' 1 , the relative peak heights (Q(3)/ 
Q(l)> Q( 5)/Q (3 ) ) for the 100 cm" 1 Stokes width compared to the 250 cm 1 
width differs by a factor of 1.37 and 0.65; i.e. the Q(l) and Q(5) lines 
are reduced in intensity by factors of 1.37 and 1.53 for 100 cm 1 width, 
thus the ratios differ by the factors 1.37 and 1/1.53. 
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III. H„ LLnmidth 

The calculation of the Raman linewidth of Hg is extremely important 
since the power of the CARS signal is inversely proportional to its square 
or for broadband CARS it is inversely proportional to the linewidth. 
Furthermore, to extract temperature measurements from the ratio of the 
signal strengths of the individual rotational lines, it is necessary to 
know the variation of the linewidth with rotational quantum number J. Un- 
fortunately, the behavior of the H 2 linewidth is complex and its variation 
with J and with temperature at high densities is not well known. 

The variation of hydrogen's linewidth with density can be divided in- 
to three regimes. At low density {<0.1 amagat) line broadening is due 
principly to the Doppler effect. For Raman scattering this line has a 
Gaussian profile whose width (FVJHM) can be written [18] 

Av D = \ Tff) [4(v| + v|v R ) sin 2 {|-) + v|] (17) 


where v 0 is the exciting frequency, v R is the frequency of the Raman res- 
onance, e is the scattering angle and the other symbols have their usual 
meaning. For CARS, only scattering in the forward direction is considered 
(8 = 0) and 


2v r ky % 

Av n = -A (2ln2 
D c m 


(18) 


In this range, the linewidth is independent of density and varies as the 
square root of temperature. 

When the mean free path for elastic collisions is comparable to the 
wavelength of the scattered light, but is less than the mean free path for 
collisional line broadening, the linewidth can become less than that due 
to Doppler broadening. This phenomena, known as Dicke narrowing [19], is 
important in reducing the linewidth of molecular hydrogen in the density 
regime between ^0.5 amagat and ^ 25 amagat [20]. Above this density, the 
lineshape becomes Lorentzian and its width increases linearly with the 
density (homogeneous or pressure broadening regime). Theoretically this 
linear variation with density is a consequence of assuming that only 
binary collisions take place. If the time between collisions is long 
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compared to the collision time, this criterion is met. Hunt et al . [21] 
have estimated for that the time between collisions is approximately 
three times the collision time for pressures of ^ 100 atmospheres for a 
hard sphere collision model. This would seem to indicate that linear 
pressure broadening would be suspect above this value. However, All in 
et al . [22] have found linear variation holds between 30 and 600 atmos- 
pheres of pressure, presumably at room temperature. Since the collison 
frequency scales n/T and the pressure as nT, where n is the number den- 
sity and T the absolute temperature, the collision frequency will decrease 
with increasing temperature for a given pressure indicating that the line 
will broaden linearly in this range at elevated temperatures. 

The effects of Doppler and pressure broadening can be accounted for 
simultaneously using an expression for the resonant part of the third 
order susceptibility developed by Henesian and Byer [23] 


x( 3 ) 


NjiicViM 

hn 2 G) 2 Acop 


I 

j 



W*(6. + iE ) 

j J J 


(19) 


where w(z) = complex error function (see Appendix E) 

= [ftj " (co 1 -oj 2 )]/ Ao> d 



Au d (= 2 ttAvq) is defined by Eq. (18), i\ is the linewidth (FWHM), and the 
other quantities are defined as in Eq. (11). Dicke narrowing can be in- 
cluded in a simple way introduced by Galatry [24]. This requires that 
l\ be replaced by 

r j - + a j( t > p < 20) 


where D 0 is the diffusion coefficient for the scattering molecule, p is 
the density, Aj is the pressure broadening coefficient, and k is an ef- 
fective wave vector for Raman scattering i.e. [20] 

k = [2k 1 k 2 (l-cose) 2 + (ki-k 2 ) 2 ]^ (21) 
which for the case of CARS (e = 0) becomes the difference between the wave 
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vectors of the incident and emitted photons. (This means that the line- 
width will be a minimum in the forward direction and a maximum in the 
backward direction. See, for example Lallemand et al . [25].) It should 
be noted that at high pressures (E.>>1) 

J 

w(z) -> — — 
i/F z 

and Eq. (19) reduces to Eq. (11), the usual expression for the suscepti- 
bility. 

The diffusion coefficient D is usually written as D 0 /p and has a 

temperature dependence that varies at T a where a varies between 1.5 

(rigid elastic spheres model) and 2 (Maxwellian models) [26]. The more 

conservative value of 1.5 is adopted here which, after the temperature 

dependence of p is removed, implies that D 0 ^T 2 . Several spectroscopic 

determinations of D 0 have been made ([20], [27]) but the measurement of 

the diffusion of parahydrogen through normal hydrogen made by Harteck and 

Schmidt [28] is the most accurate to date. Their value of 1.396 + 

2 -1 

0.0025 cm sec at 785 torr and 20° C is used in these calculations. 

There is little data on the temperature dependence of the broadening 
coefficient at elevated temperatures. The half width of the Raman lines 
is given by nvo' where v is the average velocity of the colliding pair 
and a' is the real part of the optical cross section [29]. If o' is 
assumed temperature independent, the linewidth would have a /T dependence. 
This is consistent with the data of Lallemand and Simova [30] who meas- 
ured the temperature dependence of the broadening coefficient of the Q(l) 
line of between room temperature and 150° C. (These measurements were 
made using a backward wave Raman amplifier and hence their linewidth is 
larger than that of the forward wave.) Lallemand and Simova found their 
results in disagreement with the theory of Ref. [29]. For simplicity, 
the /T" dependence of Aj will be used in our estimates of linewidth. 

There have been a variety of measurements made of the pressure 
broadening coefficients of forward direction Raman scattering of hydrogen 
at room temperature for the lowest values of 0(0-4). Table 1 compares 
these and also lists some theoretical values. Columns 1 and 2 show the 
broadening coefficients measured by spontaneous Raman scattering over a 
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TABLE 1 


H 2 Pressure Broadening Coefficients (10 3 cm* 1 /amagat) 


Line 

All in et al? 

[22] 

(1967) 

Murray 

et al. b [20] 
(1972) 

Foltz et al . c 

[34] 

(1966) 

Hunt et al . d 

[21] 

(1970) 

0wyoung e 

[27] 

(1978) 

f 

Van Kranendonk 

[29] 

(1963) 

(Theory) 

Hunt 

9 

et al . 
(1970) 
(Theory) 

Q(0) 

2.32 + 0.04 

2.2 + 0.15 


2.78 + 0.10 


2.2 

2.59 

Q(l) 

1 .40 + 0.03 

1.5 + 0.1 

1.77 

1 .83 + 0.10 

1.71 + 0.02 

0.6 

.72 

Q(2) 

2.53 + 0.03 

2.6 + 0.2 


3.05 + 0.10 


2.6 

2.87 

Q(3) 

3.66 + 0.04 

5.5 + 1 


4.60 + 0.10 


3.6 

3.87 

Q(4) 




3.26 + 0.10 


3.0 

3.86 


a. Spontaneous Raman scattering, 300-600 amagat. 

b. Spontaneous Raman scattering, 5-100 amagat. 

c. Electric-field-induced measurements, 27-55 amagat. 

d. Electric-field-induced measurements, 4.4-35 amagat. 

e. High resolution CW stimulated Raman spectroscopy, 0.26-14 amagat. 

f. Theory, isotropic scattering due to quadrupolar interaction. 

g. Calculation using theory of Van Kranendonk with revised molecular constants. 


very wide range of densities (5-600 amagat). These show good agreement 
except for the Q (3 ) lines. As noted above, Ref. [22] found the broadening 
to be linear to 600 amagat. Columns three and four display the linewidths 
measured by electric field induced dipole transitions over a more re- 
stricted density range. Brannon et al . [31] has shown that the theory of 
Raman line broadening may be applied without modification to electric 
field induced lines. These values tend to be somewhat higher than those 
measured by Raman scattering. Column 5 lists a measurement using the high 
resolution technique of CW stimulated Raman spectroscopy. Columns 6 and 7 
give theoretical values obtained using the theory of Van Kranendonk [29]. 
Hunt et al . [21] has used this theory to recalculate the linewidths util- 
izing pressure modified effective rotational constants. The disagreement 
between theory and experiment for the Q(l) line is discussed in Ref. [22], 

The pressure broadening coefficients of Allin et al . are adopted for 

the Q(0) - Q(3) lines in this study. The measured value of Hunt et al . 

is used for the Q(4) line. Finally it is useful to have a value of Aj 

for the Q( 5) line to obtain temperatures from the ratio of its strengths 

at higher temperatures. Unfortunately, there is no measured value. Using 

parameters calculated from Hunt et al., based on Van Kranendonk' s theory, 

-3 -1 

a calculated value of 1.08 x 10 cm /amagat @ 300 K was obtained. This 
is a particularly risky procedure, however, as the agreement of the modi- 
fied calculated linewidth with the experimental values becomes poor with 
higher J numbers [21]. 

Figure 3 shows a log-log plot of the Q(l) linewidth of H^ for room 
temperature and two temperatures that bracket the minimum and maximum 
temperatures that could be encountered in the SSME preburner. The density 
range is between 35 amagat (3000 psi, 1700 K) and 150 amagat (6000 psi, 

850 K) and the line broadening is well within the pressure broadening 
regime. 

As a final remark, it should be noted that St. Peters [32] has found 
that when high power lasers are used, the appropriate linewidth for CARS 
may be larger than the Raman linewidth due to the laser mode structure. 
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RAM 



Figure 3. Variation of Hydrogen Q (1) linewidth as a function of 
density for several temperatures. 
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IV. Vhah matching 

According to Eq. (7) the input laser beams should be phasematched 
in order to generate appreciable signal levels. The phasematching con- 
dition 

Ak = 2k x - k 2 - k 3 (22) 

where |k| = ~ is a vector relationship that must be simultaneously 
satisfied with the energy conservation condition 


2wi - w 2 “^3 - 0 . (23) 

If the medium has little dispersion (n x = n 2 = n 3 ) as is the case with 
gases at moderate pressures, these conditions can be met with all beams 
parallel to one another (See Figure (4a)). In this geometry the inten- 
sity of the CARS signal varies as the sine squared i.e. 


I * 




and reaches a maximum in a coherence length L where L = n/Ak. 


This 


length can be calculated by expanding the refractive index of the medi- 
um with respect to frequency [35] 

-1 


ttC 


n, an , a 2 n 
2 + 


(24) 


For gases at atmospheric pressure this distance can be tens of centi- 
meters while for liquids it can be fractions of millimeters. Since the 
CARS intensity is proportional to the square of l_ c , the signal levels 
will be unacceptably small if one works with parallel beams when the 
medium has appreciable dispersion. Hence, in these media, phasematching 
is usually achieved by crossing the input beams at a small angle (see 
Figure 4b). The interaction length in this situation is then limited by 
beam walkoff. 

It is possible to estimate the coherence length of a gas as a func- 
tion of density by calculating the frequency derivatives of the dispersion 
equation. The variation of the index of refraction with denstiy can be 
represented with a high degree of accuracy using the Lorenz-Lorentz re- 
lation [36,37] 


n 2 -l _ 


n 2 +2 


= R(x)p 


(25) 
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Figure 4. Wavevector diagrams illustrating different phasematching 
configurations, (a) Col linear phasematching, (b) Crossed beam 
phasematching necessary for media with appreciable dispersion. 

(c) Relationship between wavevectors and Ak for phase mismatched 
beams . 
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where n is the index of refraction and R(x) is the specific refraction 
that depends on the polarizability of the species. The dispersion equa- 
tion for H 2 at a density of one amagat is [38] 

no-1 = a + (26) 

where the empirical constants, a, b and X Q have been measured by Koch [39]. 
Solving Eq. (25) 


Putting p = 1 in this equation and equating it to the right hand side of 
Eq. (26) one finds R = A/(A+3) and 


where 


n" - 


3pA(x) 

= 3- ( p-1 ) A(x) 


A(x) = (a + x 2 ^z'). 


(27) 


Equation (27) is the density dependent dispersion equation. This equation 
can be simplified by noting that even at high densities (p-l)A(x)«3 (at 
200 amagat pA = 0.056 at X = 532 nm) and 

n 2 - 1 = pA(x) (28) 

Table 2 delineates the calculation of L c for the H 2 Q(l) line assuming 
a pump laser at 532 nm. This analysis indicates that the dispersion of H 2 
scales linearly with density and not with the square root of pressure as 
stated by Wilke and Schmidt [40] and except for extremely high densities 
the coherence length varies inversely with density. Figure 5 shows the co- 
herence length as a function of density. At p = 100 amagat, L c = 2 mm. 
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TABLE 2 


L 


c 


TtC 

[2n l +w 1 n"] 


where , 

n = [Ap + 1J 5 


, _ dn _ A 
dw 2n 


n" - 1 [A" 

n " n [T 

A(w) = a + b[d/u 2 -x|] 


TT p “ 

-1 



A' 


2bd 


A" 


2bd r 4d 3 ■, 

-?=■ ( ^ 2 n T ' r~p ) 


[ ] = dtjj -2 -Xo 


d = (2ttc) 2 

For the Q(l) line of hydrogen at x = 0.532 pm 

a = 2.716 (10 -4 ) 
b = 2.112 (10" 6 ) 

X§ = 7.760 (10" 3 ) for X in ym [Ref. 39] 
and 


L 


c 


1_ 1 54 

6.92(10“ 2 )f + 1.75(10“ 7 )p 2 


(cm) 


22.2 cm 

P 


p in amagat 
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COHERENCE LENGTH (mm) 



DENSITY (amagat) 

Figure 5. Hydrogen Q(l) coherence length for beams in collinear 
phasematching geometry as a function of density. 
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To predict the behavior of the CARS signal in the medium, this value 
of coherence length should be compared to the distance, b, near the beam 
waist over which a focused beam remains essentially parallel. If a beam 
of divergence half angle 6, is focused by a lens of focal length f, this 

-i 

distance is approximately given by 27tf 2 6 2 /x. Even for a conservative 

divergence of 6, =0.1 mrad, b = 12 cm when a beam passes through a 1 m 

-2 

focal length lens. Hence, at a density of 100 amagat, the CARS signal 
would have a squared sinusoidal variation with a period of ^ 4 mm in this 
region. (It is interesting to note that when b<<l_ c , collinearly phase- 
matched beams generate only 75% of the CARS signal in a length of ^ 7.5b 
[35].) 

It is clear, at the high densities expected in the preburner, the 
pump and Stokes laser beams must be crossed at a small angle. Using the 
linear momentum conservation equations implied by Figure 4b and the energy 
conservation relationship Eq. (23), an expression for crossing angle e can 
be found: 


cos e 


4o)io)2n§ - 4to£ (n <3 - nf) - (n 2 - n|) 

4o) ^ 0)2 n i n 2 


(29) 


This crossing angle was calculated for the hydrogen Q(l) line and a pump 
laser wavelength of 532 nm. The result, shown as a function of density, 
is given in Figure 6. At a density of 100 amagat the phasematching angle 
is * 0.78°. 

It is interesting to calculate the sensitivity of the CARS intensity 
to a small phase mismatch due to a non-optimum matching angle e. From 
Eq. (7) the intensity varies as sinc(x)(=sin(x)/x) , where x = Mp Ak 
is the magnitude of the wave vector mismatch (see Figure 4c), and l is 
the interaction length. Holding <j> constant at the phasematched value, 
the CARS intensity was calculated as a function of e for the Q(l), Q(3), 
and Q(5) lines of W^. The results are shown in Figure 7 for p = 1 

and p = 100 amagat density and l = 1 cm. The FWHM's of the sine 

-4 

functions are ^ 0.01 mrad (~ 7 x 10 degrees) and the intensities for 
the Q lines overlap for the low density case. However, at 100 amagat the 
increased dispersion separates the lines in 0-space. This would indicate 
that separate angles are required for the individual lines. In practice, 
however, the beam divergence of the laser beam are greater than the 
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CROSS 


PHASEMATCHING ANGLE 



DENSITY (amagat) 

Figure 6. Crossing angle required for Q(l) line phasematched beams 
as a function of hydrogen density. 


INTENSITY (ARBITRARY UNITS) INTENSITY (ARBITRARY UNITS) 



CROSSING ANGLE (DEGREES) 



CROSSING ANGLE (DEGREES) 


Figure 7. CARS normalized intensity variation for the first three 
strong Q branch lines of hydrogen as a function of crossing angle, 
(a) 1 amagat density, (b) 100 amagat density. 
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angular separation between the Q(l) and Q ( 5 ) lines and signal can be pro- 
duced for these lines simultaneously - albeit at considerably lower in- 
tensities than those calculated for the ideal case of exact phasematching. 

The details of the intensity of various lines will depend on the in- 
tensity distribution as a function of angle of the laser. By its nature 
this variation cannot be folded into the signal and error calculations of 
Section VI. Indeed, ratios of the intensities of the Q lines may depend 
on the exact angle of crossing and could differ for identical conditions 
of temperature and pressure. This is an area where experimental investi- 
gation is warranted before attempting a measurement. 


28 



V. Loaqa Intensity UjnitationA 


There are several effects that will place upper limits on the inten- 
sity of the pump beam in the medium. These in turn will limit the CARS 
signal strength and hence, the accuracy with which temperature measure- 
ments can be made. Although these effects generally are not important 
when performing diagnostics in flames at atmospheric pressure they can 
become important considerations at higher pressure. 

Two of these effects were first discussed by Reginer et al . [41]. 

The first is most important when using a col linear phasematching geometry 
and results in an augmentation of the Stokes intensity through the process 
of stimulated Raman gain. A high intensity laser beam with angular fre- 
quency oil travelling through a Raman active medium can induce coherent 
radiation to propagate parallel to it at frequency m 2 and intensity I 2 . 
Assuming negligible depletion of the pump beam intensity I x 


I 2 U) = I 2 (0) e g£l1 


[30] 


where l is the length that Ii interacts with the medium, I 2 (0) is the 
Stokes field originally present and g is the Raman gain of the medium given 
by [see for example. Ref. 42] 


(4ttc) 2 N 0 pA j / d \ 

n^JiwiT . id ft J . 

<J ' ' J 


[31] 


where N 0 is the number density at 1 amagat (2.69 x 10 -19 cnr 3 ) and all 
other symbols in this expression have been previously defined. Where the 
linewidth decreases or remains constant as the density increases, the 
gain increases (Dicke narrowing regime). In the pressure broadening regime, 
the gain is independent of density. For the Q(l) line at 1200 K and 100 
amagat, g ^ 6 x 10~ 10 cm/watt. If Ii max , the highest allowable pump beam 
intensity is considered to be that which will induce a 10% growth in the 
Stokes field (I /I (0) =1.1) then I 1 max * 1.7 x 10 8 watts cm~ 2 /£. For 
an interaction length approximately equal to the Rayleigh length (> 10 cm), 
this maximum value is quite low and presents a severe constraint that must 
be avoided if reasonable signal levels are to be attained. This points 
out another difficulty with performing the measurement using the more geo- 
metrically simple collinear phasematching. If the pump and Stokes beams 
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are crossed, this hazard can be avoided. 

It should be noted that the pump beam alone can generate a co-propa- 
gating Stokes wave with an appreciable intensity. This could yield spuri- 
ous temperature measurements even when the pump and Stokes beams are 
crossed if the CARS signal produced parallel to the pump beam were not 
separated from the crossed beam CARS signal. In this case, Eqs. (30) and 
(31) are again used to calculate the power of the stimulated Stokes wave, 
but in the absence of an external Stokes beam the initial Stokes power is 
taken as the initial spontaneous Raman noise level, i.e. P g (0) = Jiu) 2 Ij 
[42]. For the Q(l) line at 100 amagat density, this is ^ 2 x 10 -9 watt. 
Therefore, for a pump intensity of 10 9 watts/cm 2 , a distance of ^ 46 cm 
would be required to generate the same spectral power density as that 
available from the Stokes beam (assuming a 1 megawatt pulse and a 100 cm -1 
spectral width). Since the preburner diameter is 25 cm, the stimulated 
Stokes power is less than 10 -5 of this value and no special precautions 
need be considered. 

In reality, both of these calculations will overestimate the power 
generated by stimulated Raman scattering since it has been found experi- 
mentally that saturation of the Stokes intensity in begins near a 
density of 10 amagat (at room temperature) and even suffers decreases in 
intensity [43,44]. This diminution is attributed to a decrease in the 
gain at higher densities and the fact that the increased optical dispersion 
at these densities makes the four wave mixing process more difficult. 

A second phenomena that can cause measurement inaccuracies if I]_ and 
I 2 are too intense is a population perturbation between the lower and up- 
per rotational levels. Stimulated Raman coupling between these two levels 
can affect a. (see Eq. (11)) such that 

sJ 

Aj (t) = Aj(0) exp(-t/r A ) 

where, on resonance, the time constant for this perturbation is given by 



where g. = 2v. + 1 and I 2 is the Stokes intensity within the Raman line- 
width. In the density region where the Raman linewidth is homogeneously 
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(pressure) broadened, this time constant varies linearly with density. 

r. 

However for broadband CARS I 2 = I 2 tot (^-) where Aw 2 is the spectral 

width of the Stokes beam the time constant is independent of density. 

For the Q(l) of with Aw 2 = 2nc(100 cm'^ ), t a "^= (2 x 1 ) I x I 2 tot . 

If i 1 max , the maximum allowable intensity, is defined as that which causes 

a 10% population perturbation, it is assumed I 2 = 0.5I l5 and the time 

-8 

of the perturbation is taken to be the laser pulse width (10 sec) then 
Ii max „ iqIO watt CIt f^ This va i ue should not be exceeded. 

The third effect that may limit the CARS signal strength is the 
maximum pump beam power than can be focused into the H ^ without causing an 
ionization breakdown. The breakdown threshold intensity for most gases 
exhibits a minimum in the pressure range of 100 - 1000 atm [45]. The 
variation with pressure of the threshold minima is slower for molecular 
gases than for atomic gases [46]. 

Measurements of the breakdown threshold intensity have been taken for 
various gases (He, N 2 , Ar) over a large enough pressure range to determine 
the minima [47]. However, data for in the range 0.5 < P < 100 atm 
[48], do not cover sufficiently high pressures to show the breakdown mini- 
mum nor to cover the pressure region of interest to this study. The pres- 
sure at which the breakdown minimum occurs can be calculated from the 
classical impact ionization description of breakdown [45]. An approxi- 
mation to the calculated minimum should give a reasonable approximation 
for pressures near 400 atmospheres. The electron energy change is a 
maximum when the momentum- transfer collision frequency v m is equal to the 
incident field frequency. In terms of the reduced pressure (torr), 
collision probability P , and mean energy (eV), the collision frequency 
is expressed as 

v m = 5. 93 (1 0 7 ) P c \/Tpo, (p Q in torr, e in eV). 

From tabulations [49] of P c vs y/z , the mean value of P for Hg can 
be estimated to be approximately 80. It has been noted [45] that for 
hydrogen the collision frequency is very weakly dependent on electron 
energy, so 

v m - 5.9 (10 9 ) p 0 (p 0 in Torr). 
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With these two relations a bound can be placed on the pressure of the 
minimum for the data presented in Ref. 48. Figure 8a shows the ruby 
laser H 9 breakdown data (converted to power density using a beam focal 

_o 

diameter of 1.2 x 10 cm) as a function of pressure for a wavelength of 
594.3 nm (cj = 2.71 x 10^ 3 sec""'). The curve is an approximate fit ex- 
trapolated to the calculated minimum: 600 £ p Q 1 750 atmospheres. 

For this data to be applicable to this study it must be scaled to 
the expected experimental conditions. The effect of the differences in 
temperature, beam focal volume, and beam frequency are the most signifi- 
cant. Since the minimum breakdown threshold scales with density (reduced 
pressure), the data from Fig. 8a should be shown as a function of density 
(assuming the data was taken at or near room temperature). The dependence 

of the breakdown field strength on beam focal volume has been shown [45] 
-3/4 -3/2 

to vary as E « A (or S « a ' for power density), where A is the 
characteristic diffusion length of the focal volume. For a cylindrical 
focal volume with diameter d, length L, 



For most beam geometries d<<L so A « d/4.8. As the focal spot size in- 
creases, the breakdown threshold decreases. For the data of Ref. 48: 

f = 2 cm, d « 1.2 (10' 3 ) cm; so A - 2.3 (10”^) cm. Possible values 
for this study include: 

f = 100 cm, .01 <_ d <_ .05 cm; thus 2.08(10~ 3 ) £ A <_ 1 ,04(10 -3 ) cm. 

The power density focal spot scaling factor, (a/A 0 ) -3 / 2 , should be in the 
range 3.43(10“ 3 ) < (^ ) -3 / 2 < 3.73(10"^). The beam frequency dependence of 
the threshold power density should be S » u 2 , since the rate of increase 
of the electron energy is inversely proportional to the square of the in- 
cident field frequency. This naive view is consistent with the behavior 
of threshold in the frequency range of a ruby laser (used in Ref. 48) 
and a frequency-doubled Nd:YAG laser (present study) [45] 
i.e. 2.71(10^ 5 ) <_u <_ 3.54(10 15 ) sec -1 . The frequency scaling factor 
should increase the breakdown threshold by the factor (w/u) 0 ) 2 = 1.7. 

Figure 8b shows the possible range of the power density breakdown 
threshold for the experimental conditions of this study as a function of 
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S (watts cm" 2 ) S (watts cm -2 ) 




Figure 8. H 2 breakdown threshold power density as a function of 
pressure and density (a) Data from reference [48]. Curve is 
approximate fit extrapolated to predicted minimum at ^ 675 
atmospheres, (b) Expected range of breakdown power density scaled to 
parameters of interest in this study for a range of focal spot diam- 
eters. Calculated from data given in (a). 


33 






gas density. For the expected number density of 100 amagat, the limiting 

o q p 

power density will be 8.5(10 ) <_ 9.5(10 ) watt/cm , for the range of 

focal spot sizes .05 _> d >_ .01 cm. Note that similar predictions of the 
minimum threshold for single mode lasers have been within a factor of 2 in 
the power density and pressure [45,47]. If the laser is significantly 
multi -mode, the threshold will be lower than that calculated [50]. 
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VI . Signal Stating th and Expe.cXe.cl Eaaoa 


An estimate of the expected signal strength can be made using Eq. (8) 
which assumes the condition of exact phasematching (Ak=0), 


■ 


l x CARS I 


2 

II Pc 


where l is the interaction length, it has been noted that n^n^^^l and 
powers have been substituted for intensities I 2 and I 3 by assuming that 
the CARS signal is generated with the same cross-sectional area as the 
Stokes beam. The CARS susceptibility is given by Eq. (11) which on res- 
onance is 


X CARS = jfif 


2a,Nc 4 . . x , 

j / do_\ 1_ 

Vd nL r. 


where the full width r. has replaced the half width y. and it is taken 

*3 3 

that \ - 0. Combining these expressions one may calculate the number of 

NR ry, 

photons generated per laser pulse N. for a given resonance j in terms of 

<3 

the energy of the Stokes beam 


Njd.p) = 



(32) 


19 -3 

where N 0 is the number density of 1 amagat (2.69 x 10 cm ), p is the 
density in amagat, A^ is the fractional population difference between up- 
per and lower states (see Appendix C), and (^)j is the Raman cross sec- 
tion (see Appendix B). E 2 is the Stokes energy within the Raman linewidth 
which for a broadband beam of spectral width Aw 2 is 


C _ r total 
t 2 - t 2 



Hence for broadband CARS, the signal strength is inversely 'proportional to 

r. and in the pressure broadened region increases linearly with density. 

J 

The signal strength varies quadratically with pump laser intensity, an 
attractive parameter to use to increase measurement accuracy. However, 
the upper limits to the intensity set by population perturbations, stimu- 
lated Raman gain, and laser induced gas breakdown as described in Section 
V must be observed. 
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For the purposes of this calculation, a pump beam intensity of 10 

_2 

watts cm and a Stokes beam energy of 10 mJ in a spectral width of 
100 cm -1 will be considered. Modern Nd:YAG lasers are capable of pro- 
ducing nearly 100 times this intensity at a wavelength of 532 nm when 
focused by a 1 m lens. An interaction length of 1 cm is used in accor- 
dance with the discussion in Section VII. Figure 9 shows the expected 
CARS signal strength (photons per laser pulse) calculated from Eq. (32) 
as a function of temperature for the case of a constant number density of 
100 amagat for the first six Q lines of hydrogen. As is evident from 
these curves, the expected signals are quite large. The odd valued Q 
lines have greater signal strengths due to the statistical weight 
factor and hence the analysis will concentrate on these lines. 

As discussed in Section VII, there are several methods that can be 
used to recover and analyse the CARS signal. Two are considered here. 

The first is to ratio the signals obtained from the individual lines and 
the second is to ratio the signals obtained from two interference filters. 
Figure 10a shows calculated ratios of the odd (strong) lines using 
Eq. (32) as a function of temperature for a constant density of 100 amagat. 
Figure 10b shows the expected ratios if the CAR spectra is convoluted by 
1.0 nm (FWHM) interference filters with triangular spectral profiles 
(characteristic of two period filters). The filters are centered at 438, 
437, and 436 nm respectively (see Figure 19). Note that these two figures 
have different ordinates. 


The Q 5 /Q-| ratio has the greatest slope of the Q ratios and varies more 
than an order of magnitude for the range of temperatures of interest 
(800-1600 K). The F^F^ ratio has the greatest ratio among the F ratios, 
but the slope is considerably less than that of the Q5/Q-] » its value 
varying less than a factor of 4 in this same temperature interval. 


Error Analysis 


It is useful to consider the accuracy with which the temperature can 


be extracted from these ratios. This can be done in terms of simple error 


propagation theory [51]. The variance a l of a quantity x in terms of the 

A 

variances of variables u and v is approximately 


a 


2 

x 
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Figure 9. Calculated CARS signal strength for Q branch lines of H 2 at 100 amagat density. 
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Figure 10. (a) Calculated CARS Q line signal ratios for 100 amagat 

density, (b) Calculated ratios from CARS H 2 Q line signals taken 
through the three interference filters displayed in Figure 19. 






where the covariance o 2 v is zero if the fluxuations in u and v are un- 
correlated. The error in temperature due to error in the ratio of counts 
R is simply 


«f ■ 


(£L> 


R ' 3R 


where R = N!/N 2 . The uncertainty in the ratio will consist of the experi- 
mental uncertainty with which the ratio can be measured a 2 R (due to 
gated integration electronics, ratiometers, A/D etc.) and that due to 
uncertainty in the number of counts a so 

NR 

a 2 = a 2 + a 2 
R MR NR 


and 


,3Rn 2 


,2 = a 2 f £!!)“+ a 2 (3R.1 
? NR °Ni 3Ni a N 2 ' 3N 2 


The uncertainty in the number of counts can be ascribed to purely statis- 
tical causes a ^ ^ and to the dependence of this number on the linewidth, 
which itself has appreciable uncertainty. If the measurements are taken 
using the interference filter configuration, the counts will depend on 
contributions from several lines so 


a 2 = a 2 + 
N stat 


V n-2 v 2 


Then, noting ]||| = ^ and letting a = \ZTT the variance in temperature 
due to these uncertainties can be written 


a 


2 

T 




N 2 + Ea r 2 



Figures 11 through 16 show the calculated error aj for the same con- 
ditions used for the signal strength plot (Figure 9). Figures 11 and 12 
consider uncertainty in temperature due to statistics, ratio measurement 
and linewidth uncertainty (two cases) individually. These calculations as- 
sume that only 5% of the photons are recovered and impinge on a photomul- 
tiplier with a quantum efficiency of 0.2. For these conditions, Figure 11a 
shows that the projected error in temperature due to statistical uncertain- 
ty alone is quite small. Figure lib displays the uncertainty in T expected 
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Figure 11. Calculated temperature standard deviation for Q line ratios, 
(a) a, due to signal statistics only, (b) due to 0.5% projected error 
in measured signal ratio. 





due to measurement error in the ratio. These curves were calculated as- 
uming a 0.5% relative accuracy (a /r = 0.005), a reasonable value for 
matched detectors when the variation of R is less than two orders of 
magnitude. 

The uncertainty in the linewidth broadening coefficient can lead to 
appreciable error in the temperature at high densities. Figure 12 dis- 
plays the error for two cases. If the listed uncertainties of Allin et 
al. [22] are used for lines Q(0) - Q(3) (column 1 Table 1) and 0.1 x 10~ 3 
cm“Vamagat is the assumed error for lines Q(4) and Q(5), the resulting 
error in T due to this linewidth uncertainty alone is shown in Figure 12a. 
However, the spread among the broadening coefficient values listed in 
Table 1 is considerably greater than this error. Figure 12b displays 

a T due to the uncertainty in r. above if the various broadening coeffi- 

1 -3-1 

cients are assumed to have constant uncertainty of 0.25 x 10 cm /ama- 

gat. Clearly, this uncertainty in the linewidths will contribute sub- 
stantially to the uncertainty in the temperature if the published values 
in the linewidth are used in the analysis of the data. Using Eq. (31) the 
total calculated uncertainty in the temperature can be computed. The re- 
sults, using the above parameters and the smaller uncertainties in the 
line broadening coefficients, are shown in Figure 13. 

It is interesting to note that the ratio that yields the least uncer- 
tainty differs for each of the plots. Under the assumptions of Figure 13 
the preferred ratio is the Q 3 /Q-j in accordance with the linewidth uncer- 
tainties used in Figure 12a. This would not be the case for the more 
liberal uncertainties assumed for Figure 12b which indicates the Qg/Q-j 
ratio is the preferred quantity. Near the operating temperature of the 
preburner (^ 1200 K), these calculations indicate a temperature uncertain- 
ty of ^ 10-40 K (depending on the linewidth uncertainty assumed) or an 
accuracy of 1-3%. 

It is possible to extract temperature from the measurements of the 
ratios without using the linewidth data by precalibrating the instrument 
using a variable temperature high pressure cell. The uncertainty would 
then approach that given by Figure lib. That is, it would be due to the 
experimental uncertainty in the ratio. In this case, the Qg/Q-j ratio 
would be the preferred value to use for temperature extraction consistent 
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Figure 13. Calculated temperature standard deviation for Q line ratios when the effects 
of statistics, the linewidth uncertainties of Ref [22], and ratio measurement uncertainty 
are included. 



with the fact it has the steepest slope. Then the temperature uncertainty 
would be ^ 1.5 (K) at 1200 K corresponding to a 0.1% accuracy. 

Figures 14, 15 and 16 show the corresponding error calculations for 
the filter ratios. Here, the calculation assumes that 40% of the signal 
is recovered and is divided by a 50% beam splitter, passes through an 
interference filter with 50% peak transmission and is focused onto a PMT 
with 0.2 quantum efficiency. The results are similar to those of the Q 
ratios. Again the uncertainty due to statistics is negligible and some- 
what less than those of the Q ratio due to the higher signal throughput and 
the fact that the filters accept contributions from several lines. The 
linewidth uncertainties yield somewhat higher temperature uncertainties 
for the filters because the uncertainties from several lines contribute 
to the signal from each filter. Using the ratio. Figure 16 projects 

a total uncertainty of ^ 25 K, about 2.5 times greater than that for the 
Q ratio. Considering the case of a calibrated instrument and using the 
uncertainty implied by the measured ratio accuracy of Figure 14b, the 
R^/R^ ratio projects only a slightly greater uncertainty in temperature 
(1.75 K) than the Qg/Q-j ratio. 

It should be noted that there is another source of error that could 
not be included in these calculations. A variation of density in the pre- 
burner due to temperature and/or pressure variations can change the phase- 
matching angle. Not only will this have the effect of reducing the signal 
strength, but could also change the relative intensity of the spectral 
lines in a random fashion as discussed in Section IV. Furthermore, the 
CARS signal strength used in these error calculations consider the case, 
of perfect phasematching. This cannot be the accomplished for all the Q 
lines simultaneously and the signal levels obtained here are surely over- 
estimates. This would imply greater statistical errors than those shown 
in Figures 11a and 14a. 
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Figure 14. Calculated temperature standard deviations for 
filter ratios, (a) due to statistics only, (b) a y due 
to a 0.5% projected error in measured signal ratio. 
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TEMPERATURE (K) 

Figure 15. Calculated temperature standard deviation for 
filter ratios, (a) a T due to linewidth uncertainties of 
Ref [22]. (b) ct t due to constant linewidth uncertainty of 

2.5 x 10“ 4 cm -1 /amagat. 
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Figure 16. Calculated temperature standard deviation for filter ratios when the effects 
of statistics, linewidth uncertainties of Ref [22], and ratio measurement uncertainty are 
included. 



VII. JnA&uumnt and Facltity Con&idtiAa£loni> 


In this section, the experimental layout of a CARS system is briefly 
discussed. This system would be capable of producing laser beams and CARS 
signals with the characteristics discussed in previous sections. 

Generally, a CARS system can be divided into three parts: a source of 
high intensity radiation, beam transfer and focusing optics, and signal 
recovery and processing assembly. The beam transfer optics are neces- 
sarily dependent on the physical design of the test facility and will be 
discussed last. Due to signal strength considerations the choice of laser 
for the pump radiation to* is limited to high intensity pulsed lasers for 
most practical purposes. Of those that are readily available, ruby and 
Nd:YAG lasers offer the highest intensity. The latter is the first choice 
since its doubled output can pump visible dyes useful for creating the 
Stokes laser radiation and has the advantage of offering higher pulse rates 
resulting in a higher data acquisition rate. 

Figure 17 displays the optical design of the front end of a CARS 
system built around a Nd:YAG laser which is suitable for ^ temperature 
measurements. A portion of the second harmonic of the Nd:YAG laser serves 
to pump a broadband dye laser oscillator-amplifier system. After passing 
through telescopes that modify the beam diameters and divergences, the 
pump and Stokes radiation are rendered parallel and relayed to the meas- 
urement point. The characteristics of this and similar systems have been 
discussed previously [52, 53]. Shirley et al . [7] have done single pulsed 
CARS thermometry of in laboratory burners using a similar set-up. They 
find a dye mixture of cresyl violet perchlorate and oxazine 725 satisfac- 
tory to efficiently produce the Stokes radiation at 681.5 nm with an 
acceptable bandwidth when pumped by the second harmonic of a Nd:YAG laser 
(2xNd). This set-up is sensitive to misalignment and must be placed in a 
mechanically and thermally stable environment. 

There are several methods to process the recovered CARS signal. The 
three considered in this study are shown schematically in Figure 18 in 
order of decreasing complexity. The first two require monochrometers 
while the last needs only narrow band interference filters. The first 
method will yield the greatest information. Here all the lines of a CAR 
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Figure 17. A suitable design for the front end, or laser sources, of a 
CARS system. Code: P, prism; 2X, KDP doubler; M, mirror; D, dichroic, 

BD, beam dump; IRF, infrared filter; DC, dye cell; L, lens; BS, beamsplitter; 
OF, optical flat. 
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Figure 18. Three variations of detecting and processing H 2 CAR 
spectra, (a) Viewing complete spectra using optical multichannel 
detector (OMD). (b) Viewing individual odd numbered Q branch 
lines, (c) Viewing portion of spectra through interference 
filters. 
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spectrum within an interval defined by the spectral width of the broad- 
band Stokes laser are viewed using an optical multichannel detector (OMD). 
This is the most common arrangement used for N 2 CARS thermometry. How- 
ever, the OMD is expensive and very sensitive to mechanical vibration. 

The second technique elimates the OMD and collects the signals from 
individual Q branch lines using fiber optics and transfers these to photo- 
multiplier tubes. A one meter spectrograph using a 2400 l/mm grating can 
produce a dispersion of % .3 mm/A at the anti-Stokes wavelength. This 
results in a spatial separation of ^ 1.8 mm between Q(l) and Q(3) lines 
and % 3.3 mm between Q(3) and Q(5) lines - a sufficient distance to allow 
the placement of optical fibers. The electronic signals are then pro- 
cessed and ratioed. 

The last method considered eliminates the spectrograph and uses two 

o 

interference filters. Interference filters with bandpasses of 10 A are 
commercially available. Figure 19 shows the spectral position of filters 
with triangular passbands (somewhat characteristic of two period filters) 
used in the calculations of Section VI with the positions of several of 
the Q branch lines indicated. As noted in that section, this technique 
is somewhat less sensitive to temperature changes than taking ratios of 
the signals from individual lines. However, due to its simplicity, 
mechanical stability and insensitivity to signal beam position, it is 
experimentally the method of choice. This technique is only marginally 
less accurate than system (b) of Figure 18 when using a ratio calibrated 
CARS system. 

The Fuel Preburner and Test Facility at NSTL 

Figure 20 shows vertical and horizontal cross sections of the SSME 
fuel preburner. A photograph of the interior of the fuel preburner is 
shown in Figure 21a. The. interior diameter of the main wall is 10.89 
inches and surrounds a thin (0.05") liner with a 10.43" interior diameter. 
On some versions, there are two ports of ^ 0.4" diameter located ”60° 
and + 1 65 from the 0° degree line pointed toward the preburner vertical 
axis (see Figure 21b). These ports may or may not have access to the in- 
terior gas flow through the liner. Three vertical baffle plates extend 
down from the injector head to the centerline of ports at 120° inter- 
vals. 
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Figure 21. Two views of SSME fuel preburner, (a) Interior view looking 
towards injector head, (b) Exterior view of ^1.0 cm diameter port. 
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At NSTL the preburner is located on the 4th floor of the test stand 
building. During assembly the SSME is surrounded by temporary metal 
flooring which allows access to the engine components. During a test, 
these grates are removed and the remaining flooring is ^ 15 feet from the 
SSME. This permanent deck, mounted on a n, 2% feet thick concrete wall, 
is one level below a room that houses the electronics for processing the 
data received from the sensors placed on the test engine. This amplifier 
room is a vibration-free, thermally stable environment. 

A possible physical configuration for the location of the laser sys- 
tem and beam transfer optics is shown in Figure 22. The CARS system 
shown in Figure 17 could be placed in the amplifier room 1 near the con- 
crete wall. The parallel pump and Stokes beams would pass through the 
floor of this room, be rendered horizontal by a prism secured to the 
concrete wall, pass through the wall and be transfered to a mirror mounted 
directly to the preburner by another mirror mounted directly to the other 
side of the concrete wall. Because of the severe vibration, no optics 
could be mounted on any other surrounding structures. 

Figure 23 shows a mechanism that could accept the transferred beams. 
Here a lens that rides on a travelling nut of a lead screw assembly serves 
to focus and cross the incoming beams. The crossing point of the beams 
defines the region where the measurement is taken and this point is trans- 
ferred across the preburner by the movement of the lens. The CARS signal 
thus created passes through a spectral filter and is focused by a short 
focal length lens onto an optical fiber. This fiber transfers the signal 
back to the amplifier room for processing. The focal length of this lens 
and the diameter of the fiber are chosen such that maximum diameter of the 
focused signal beam does not extend beyond the diameter of the fiber optic. 
The ports are shown in the same position as those that currently exist on 
the fuel preburner. The 12" to 14" clearance necessary for the travelling 
lens assembly is available in the current setup. Only ^ 3 inches of 
clearance is available outside the other port and hence the need for a 
compact recovery arm. With ports placed at these positions, a temperature 
profile would be taken across a chord instead of the diameter of the pre- 
burner. This is a desirable situation from an experimental point of view 
as the gas streamlines are expected to be static in this region while the 
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Figure 22. Fifth floor plan of test stand facility at NSTL. The 
possible location of a CARS system is shown in the amplifier room 
which is one level above the fuel preburner. Positions of mirror 
mounts attached to the concrete wall are indicated. 
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Figure 23. Hypothetical beam receiving arm attached to fuel preburner. All optics and beam paths are 
hermetically sealed. 



center is believed to be subject to turbulence. The windows of the ports 
could be made of crystal ine sapphire as it has excellent mechanical 
strength, high temperature stability, chemical inertness, and good op- 
tical transmission (its melting point is above 2300 K and a 0.5 inch diam- 
eter, 0.2 inch thick window has a hydrostatic burst point of 9600 psi with 
a built-in safety factor of 3 [54]). Because of the cryogenic fuel flowing 
through the rocket system, all fixtures attached to the preburner must be 
purged to avoid frozen condensate. Furthermore, all beams must be trans- 
ferred in hermetically sealed tubes as the air adjacent to the SSME be- 
comes clouded during tests. 

It is useful to consider a few numbers in light of the geometry of 
the fuel preburner and test facility. As discussed in Section IV, hydrogen 
has considerable optical dispersion at high densities requiring a simple 
two beam phasematching geometry. At a density of 100 amagat the crossing 
angle 0 is ^ 0.8°. Using simple geometry, two beams of width d can over- 
lap over a length of no more than l = d/sin(e/2). The diameter of a laser 
beam with a Gaussian intensity profile focused by a lens of focal length 
f is 'v 2f6 1 / 2 where &-l/ 2 is the half angle divergence of beam. As indi- 
cated in Figure 23 a lens with a 36 cm focal length would allow measure- 
ments to be made across a chord of the preburner. Taking <$i/ 2 = 0.25 mrad 
(factory specification of a Quanta Ray Nd:YAG laser) one finds l a. 2 cm. 
This crossing angle and lens also result in a pump and Stokes beam separa- 
tion of ^ 0.5 cm at the lens. If both beams have a 0.5 cm diameter, then 
the ports in the preburner should have approximately 1 cm diameter (or 
somewhat larger to facilitate alignment). In accordance with the dis- 
cussion in Section IV, the angle between the Stokes and pump beams must be 
maintained to an accuracy better than the individual beam divergences. 

_3 

This indicates that beam separation tolerance of less than 8.75 x 10 cm 
over the path length of the relay system (<\, 25 meters). Our experience 
is that this can be achieved if single optical components are used to re- 
lay both beams. However, these components must not be stressed and the 
vibration of these pieces must be less that required to cause the focused 
CARS signal to wander from the end of the fiber optic (^ 1 mm for a 
macrofiber). Furthermore, it is likely that once set up the system will 
require an initial alignment under operating density to maximize the 
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signal . 

In order for this design to be implemented, the following conditions 

must be met. Two facing ports with at least 1 cm diameters and an unob- 
structed line of sight view of each other must be available in the pre- 
burner. In current preburners, the ports have marginally adequate diam- 

eters, but do not face each other and offer no through-beam access. 
Furthermore, half the area of the current ports are obstructed by the 
baffles extending from the injector head. Line of sight access must exist 
from the main concrete wall to the preburner, through the concrete wall, 
and from the outside of the concrete wall to the amplifier room. It is 
believed that this access could be made available, albeit with some modi- 
fications to the test facility. Mirror mounts and housings must be 
anchored directly to the concrete wall and not to surrounding structures. 
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VIII. ConcZuAionb and RzcommndationA 


The theoretical analysis in this report indicates that temperature 
measurements can be made using CAR spectra of hydrogen with the spatial 
and temperature precision indicated in NASA RFP No. 8-1 -2-EC-25732 

o 

(0.25 to 0.5 in and 10 °R). The temperature resolution for a calibrated 
instrument will be limited by the precision with which the ratio of sig- 
nals from individual rotational lines can be measured. Spatial resolution 
can approach ^ 2 cm in the longitudinal direction (with respect to beam 

direction) and ^ 0.01 cm in a tranverse direction. The time resolution 

-8 

would be limited by the pulse width of the laser (^ 10" sec) with a 
data acquisition rate limited by the pulse rate available from high power 
Nd:YAG lasers (currently ^ 30 Hz). 

A study of the physical environment of the test facility and access 
to the currently configured fuel preburner does not present such a san- 
guine picture. As presently configured, there is no acceptable optical 
access to the preburner. CARS requires two line of sight optical ports 
with clear access diameters of ^ 0.5 in. and modifications to the pre- 
burner would have to be made. Furthermore, the area adjacent to the SSME 
is subject to extreme vibration during a test and is unsuitable for the 
placement of spectroscopic equipment. This requires a system of relay 
optics over which critical alignment must be maintained. 

Because there are unknowns in the physical environment during a test 
firing, several on site tests are recommended before proceeding with a 
CARS temperature measurement. The amount of physical motion experienced 
by optical mounts anchored to the main concrete support wall and an arm 
attached to the preburner port must be measured. The observation of 
beam wander over the projected light path during a test firing measured 
using a HeNe laser and a reticon detector would be a simple and useful 
test. In addition, the amount of deviation experienced by a laser beam 
passing through the preburner during a firing also should be measured. 
Significant deviations caused by turbulence-induced index of refraction 
gradients could cause significant dephasing of the pump and Stokes beams. 

Since there are several untested aspects of the CARS technique 
in this environment it is also advisable to undertake a number of 
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laboratory measurements using a variable temperature high pressure cell. 
CARS generated in high pressure hydrogen would allow a system to be cali- 
brated taking note of the sensitivity of the variation of line intensities 
as a function of the phasematching angle. These studies could also be 
used to eliminate many of the uncertainties associated with the linewidth 
pressure broadening coefficients considering the paucity of data avail- 
able on their variation as a function of temperature. In addition, this 
work could set experimental limits on the effects of stimulated Raman 
gain and laser induced gas breakdown for hydrogen. Finally, more real- 
istic limits on the absolute signal strength could be obtained since 
the theoretical limits are seldom realized even under the most ideal 
experimental conditions. 
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Appendix A 

Density to Pressure Conversion 

The conversion from pressure to density for a given temperature is 
readily accomplished using the equation of state. Using Van der Waal's 
equation for H 0 [55] . , 

2 (p ♦ J$.) (v - nb) . nRT 

where n = number of moles of gas 

R = 0.08206 liter atm mole 1 K” 1 
P = pressure in atmospheres 
V = volume in liters 

and for 

p 

a = 0.2444 (liters) atm/mole 
b = 0.02661 liters/mole 

If n e n/V (molar number density), then the above equation can be written 

abn 3 - an 2 + (bP + RT)n = P 

and n/no is the density in amagat where no is the molar number density at 
T = 273.15 K and P = 1 atm. Figures A1 and A2 show pressure verses density 
for the temperatures, pressures and densities important to this study. 
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Appendix B 
Raman Cross Section 


The Raman cross section for forward scattering as a function of 
rotational quantum number is [14] 


da 

dn 


Q 





(v + 1) 


(Bl) 


da 

dfi 



J_ _z_ 

2o)o 45 




(v + 1) 


It is often useful to be able to calculate a and y from the Raman 
scattering cross section a zz (the cross section for scattering an in- 
cident beam travelling along the x-axis with polarization in the z- 
direction into the x-y plane with polarization in the z-direction) and the 
depolarization p (the ratio of light intensity scattered in the y-direction 
polarized in the x-direction to that polarized in the z-direction for an 
incident beam travelling in the x-direction polarized in the z-direction). 


a 2 _ 2wo (c \ 4 (3 - 4p) „ . x 2 _ 30u) o (c \ 4 p „ 

M ~ 3jig y ^> 2 j 1 + p a zz ’ M \u> 2 ) 1 + p a zz 


Taken from 


p = 45a /3 + 4y 2 ; a zz = 90Mto 0 (w 2 ) 9 V ^ 45 “ 2 + 7 
See Ref. (56, 57, 58). 

The Placzek-Teller coefficients are available from several references 
(59, 60, 61). For a diatomic molecule 

kJ - J(J + D 

b J ' (2J-1 )(2J+3 ) 

U J _ 3 (J+l ) ( J+2) 
u+2 2 (2J+l )(2J+3) 

.J _ 3 J(J-l) 
d J-2 " 2 (20+1) (2J-1 ) 

g v = v + 1 is contributed by the vibrational matrix element. 
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Application to Hydrogen 

(o 2 = 2.714 x 10 15 rad/sec (14396 cm -1 ) 

to = 7.832 x 10 14 rad/sec (4155 cm" 1 ) 

0 

p = 0.07 

2j- 2 = 1.8 x 10” 8 cm 4 gm" 1 ) 

M Ref [10] 

v 2 _Q 4 _] l 

= 1 . 08 x 1 0 0 cnr gm J 

and i o o _i 

= 9.75 x 10" cm erg [17] 

NR 

Analysis shows [62] that for very light molecules the Q branch matrix 
element is subject to a J-dependent correction factor i.e. Eq.(Bl) must be 
multiplied by (see [63]) 

( 2B V 

where a 1 is the first Dunham anharmonic coefficient [64]. Using Eq . (1 9) 
of Ref [65] a 1 can be written 

H o ^ 

a — §. _ 3 

a l T6Bf 13 
e 

= -1.62293 cm" 1 

where the remaining constants are defined in Appendix D. 
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Appendix C 

Population Difference Between Initial and Final States l . 


The fractional population difference between upper and lower rota- 
tional-vibrational energy levels is 

A j = f v'J' f vJ 

where f yJ is the fraction of molecules in the energy level designated by 
v and J. 

In thermal equilibrium, the fraction of molecules in a quantum mech- 
anically allowed eigenstate is equal to the Boltzman factor divided by 
the partition function Q [66] 


f. 


i 


exp(-E.j/kT) 

Q 


If there are degenerate states, the fraction of molecules in a given energy 
level is weighted by the degree of degeneracy or the statistical weight g k 

g k exp(-E k /kT) 
f k = Q 


where the partition function is given by 


Q = Z exp(-E./kT) 

all 

states 


= Z g k exp(-E k /kT) 

all 

levels 


Considering all degrees of freedom (electronic motion and spin, nuclear 
motion and spin) but considering A-doubling to be zero and all electron 
and nuclear spin levels to be degenerate, f j can be written as 

(2 - 6_.)(2S + 1) (2J + 1)[Z(2T + 1)] 
f vJ - Q exp(-E vJ /kT) 


where (2 - 6 QV ) accounts for A-doubling 

(2S + 1) is the electronic spin statistical weight 
(2J + 1) is the statistical weight of a given rotational level 
excluding nuclear spin 

e( 2T + 1) is the sum over all allowed nuclear levels weighted by the 
number of their magnetic sublevels. 
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The values of T are taken from the individual nuclear spins of the 
atoms (I ls I 2 )» i.e. 

T = | I i + I 2 |.|Ii + I 2 - It* ••• III - I 2 I 

However, for homonuclear molecules, the symmetry of the total wavefunction 
is preserved by allowing only even values of T for symmetric wavefunctions 
and only odd values of T for antisymmetric wavefunctions. 

If the various degrees of freedom are considered to be completely de- 
coupled (which is not always the case, as indicated below) then the par- 
tition function can be approximated by 

0 = 0,0 Q J Q 
H y el y v y rot y es y ns 

For all temperatures important to combustion, only the ground elec- 
tronic state will contribute appreciably to the sum (ignoring A-doubling). 

E e l 50 so Q el = 

The vibrational term energy can be written 

G(v) = u> e (v + iy) - “ e x e (v + j) 2 + ... 

Defining the zero point energy (u /2) to be at E = 0 

G(v) = vU - u xv + u y v 2 + a) z v 3 + ...) 
e e e e e e e 

v = 0, 1 , 2, . . . 

Q v = 1 + exp(-hcca e /kT) + exp(-2hcw e /kT) + ... 

The rotational energy can be approximated by 
F(J) = BJ(J+1) - DJ 2 ( J+l ) 2 
J = 0, 1, 2, ... 

Q rot = Z exp(-Ej/kT) = Z (2J+l)exp(-F(J)/kT) 

stiles levels 

= (1 + 3 exp(-2hcB/kT) + 5exp(-6hcB/kT) + ...) 

For sufficiently large T or small B (see Ref. 67); 
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Q rot “ f^ 2J+1 ) exp[-hcBJ(J+l)/kT]dJ = ^ 

o 

Q es = I exp(-E es /kT) 
all 
states 

= (2S + 1) for degenerate electron spin states. 

Q„. - I exp(-E/kT) - I (2T + 1) exp(-E /kT) 


all 

states 


all 

levels 


= £ (2T + 1) for degenerate states. 
T 

The Boltzman factor is given by 

f B0LTZ = exp E' E T0TAL /kT J 


E T0TAL - he T e +G(v)+F(J) 


T g = 0 for ground electronic state (cm~^) 


G(v) = G(v) - G(0) 




ai Z. 


- r/..\ 1 / “e x e j_ re , "e el 

- S(v) - 2 \^e ~ — + T" + “8 

(= 0 for lowest vibrational level) 


(cm -1 ) 


F(J) = BJ(J+1 ) - DJ 2 (J+1) 2 
Application to Molecular Hydrogen 
Consider the ground state of 

xl£ g 5 * 0 
A = 0 

I = 1/2 + T = 1 , 0 

even J -> +, s (symmetric under electron and nuclear 
coordinate inversion) 

odd J a 

Then 

Q “ (1) [1 + exp(-hcui e /kT)](2s + 1) 

where the rotational-vibrational part of the homonuclear molecule cannot 
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be decoupled because of symmetry considerations, i.e. 

Q = £ exp(-E/kT) = £ (2T + 1)(2J + 1) exp(-E /kT) 

^ all allowed allowed T,J 

states 

= £ (2J + 1) exp(-hcFj/kT) + 3 £ (2J + 1) exp(-hcFj/kT) 

even J odd J 

= £ (2J + 1) exp(-hcFj/kT) + 2 £ (2J + 1) exp{-hcFj/kT) 

j even J 

since only even T values are allowed for even J states. The above series 
can be truncated when the exponentials become sufficiently small. Then 


and 


where 


Q(T) “ 


1 + exp 


7-6331 


Q (T) 

X RN 



g T (20 + 1) 

Q(T) exp 

g T (2J + l) 

Q(T) ~ 6XP 


[tT (T e + ^ (v) + F(J)) ] 

1"21 70 - G(v) - F(0)1 
. 695T 


= 1 for even 0 
= 3 for odd J 
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Appendix D 

Energy Differences fij and Line Positions 

The difference in energy between upper and lower energy levels par- 
ticipating in transitions 0 v j is 


where 


^vj E v'j' E v"j" 


E vj = hc(T e + G(v) + F s (j)) 

T g = electronic energy (cm - ^) 

= 0 for ground state 

1 \ 2 , .. . U 3 


G(v) = u g (v + Jr) - u e x e (v + \) + “ e y e (v + Jr) + w e z e ( v + \) + 


F S (J) = B y 0(0 + 1) - D V J 2 (J + l) 2 + H v J 3 (J + l) 3 


B 


v 


D 


v 


= B e - a e (v + 1 ) + Y e (v + ^ 
= D e + B e (v + 1) 


The subscript s differentiates electron spin components. For multiplet 
states there is an expression for each spin component characterized by a 
spin splitting constant. The form of the expression depends upon the 
angular momentum coupling employed. 


For H 2 (all values in cm" ^ , Ref 34) 


u = 4401.213* 
e 

u x = 121.336* 
6 6 

Ve = °- 8129 * 
Ve = °*° 


B = 60.8530 
e 

a 0 = 3.0622 + 0.0577 (v + h) + 0.0051 (v + h) z 
D y = 0.0471 - 0.00274 (v + h) + 0.00040 (v + h) 1 
H v = [4.9 - 0.5 (v + %)]10" 5 


* accurate for v = 0, 1, 2, 3. 


In addition, each of the line positions is pressure shifted. To 
good approximation this shift can be written for the Q branch lines as 

Q(J) = ajp + bjp 2 
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where p is the density and b, is small and will be taken to be zero 

U J.L 

here, aj is related to the relative population of the initial J tn level 
n,/n [68] 

n J 

a J = a ° + a c T- 

For a 0 = -(2.0 ± 0.2) x 10~ 3 , a = -(1.7 ± 0.4) x 10~ 3 cm"Vamagat at 
300 K. 
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Appendix E 

Analytical Form of Anti -Stokes Spectrum 

The formula for the general case of a broadband Stokes laser of 
Gaussian profile and narrow, but finite band pump laser beams is given 
by Eq. (55) of Yuratich [12]: 

S as (w as ) = [ ex P(-<S a s / ^d )/(ir Y 2 d) ] 


x 




a. % 
Y2 Y d 


Xx a.Im(w.) 

^ A NR ] ' J 


where w. 
J 


w 


n 2d 5 as 


A. + 
_J 


1Y, 


'Xj 'Xj .'Xj 


v w t \j 

W Y 2d 



Y 2 Y d 


X a .a. Im 
j,k k 



5 j + 5 k 


+ i ('v + 


Y. 


'V. 

Y. 


and w(z) is the complex error function given by 


w(z) = exp(-z 2 ) 


/ co 

exp(- 

z-t 




1 + 


dt, Imz > 0. 


- f 

Jo 


exp(t 2 )dt 


, = 4 . da 

j hn^f A j d a 


Nn^c 


6 


as 


“ M 


as 



Y 2d 



Y d “ (Y? + Y 2 ) 2 

n 2d = (1 + y|/yP _ 1 = (Y d /Y 2d ) 2 


A kj = \ - A "j ; A j = - K “ a V 

? - y/ Jffi 

u° , o)°, id? = the anti-Stokes, Stokes, and pump center angular 
as s i 

frequencies, respectively. 

Yi, Y2 j y s , Yj = the pump #1, pump #2, Stokes, and Raman 
linewidths (HWHM), respectively 

n. = Raman resonance angular frequency (calculated in Appendix D). 
3 
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Appendix F 

Computed CAR Spectra 

The following figures show calculated H 2 CAR spectra for a range of 
temperatures: 800 < T < 1800 K, AT = 100 K. Each figure consists of 

four spectra at the same temperature but computed with differing number 
density and Stokes beam width: p = 1 and 100 amagat; r = 100 and 

-i s 

250 cm . The other fixed parameters assumed for the computation are: 


Stokes beam center 

14679 cm -1 

(5812 A) 

pump beam center 

18797 cm -1 

(5320 A) 

pump beam width (FWHM) 

.8 cm”^ 



Each spectrum is normalized to the maximum peak intensity of that spectrum. 
Instrument profile convolution is not included. 
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F. 2. Computed H„ CAR Spectra at T = 900 K. 
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F. 3. Computed H 9 CAR Spectra at T 
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F. 4. Computed H 9 CAR Spectra at T = 1100 K. 
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F. 5. Computed CAR Spectra at T = 1200 K. 
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F. 6. Computed H 2 CAR Spectra at T = 1300 K. 
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F. 8. Computed CAR Spectra at T = 1500 K. 
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F. 9. Computed CAR Spectra at T = 1600 K. 
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